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ABSTRACT 

We study the effects of ellipticity in lens galaxies and external tidal shear from neighboring objects on the 
statistics of strong gravitational lenses. For isothermal lens galaxies normalized so that the Einstein radius is 
independent of ellipticity and shear, ellipticity reduces the lensing cross section slightly, and shear leaves it 
unchanged. Ellipticity and shear can significantly enhance the magnification bias, but only if the luminosity 
function of background sources is steep. Realistic distributions of ellipticity and shear lower the total optical 
depth by a few percent for most source luminosity functions, and increase the optical depth only for steep lu- 
minosity functions. The boost in the optical depth is noticeable (>5%) only for surveys limited to the brightest 
quasars (L/L* > 10). Ellipticity and shear broaden the distribution of lens image separations but do not affect 
the mean. Ellipticity and shear naturally increase the abundance of quadruple lenses relative to double lenses, 
especially for steep source luminosity functions, but the effect is not enough (by itself) to explain the observed 
quadruple-to-double ratio. With such small changes to the optical depth and image separation distribution, 
ellipticity and shear have a small effect on cosmological constraints from lens statistics: neglecting the two 
leads to biases of just ASl^ = 0.00 ± 0.01 and A£7a = -0.02 ± 0.01 (where the errorbars represent statistical 
uncertainties in our calculations). 

Subject headings: cosmology: theory — gravitational lensing 



1. INTRODUCTION 

A circularly symmetric gravitational lens is a useful theo- 
retical construct that, most likely, will never be observed in 
a cosmological setting. 7 Every real cosmological lens will 
have some small asymmetries either in its own mass dis- 
tribution (e.g., ellipticity), or in the distribution of objects 
near the line of sight (leading to a tidal shear). In fact, it 
is well known that ellipticity and shear cannot be ignored 
in models of individual observed strong lens system s (e.g., 
iKeeton. Kochanek. & Seliakll997HWitt & Maoll99l . 

Nevertheless, most analyses of the statistics of grav- 
itational lenses have used symmetric lenses. The sta- 
tistical calculations offer enough intrinsic challenges that 
most authors have stuck to idealized spherical lenses, such 
as the singular isothermal spher e (SIS) or the general- 
ized Navarro-Frenk- White (GNFW; Nava rro. Frenk. & White! 
19971 laiaol fl99^) profile fe.g.. i Narayan & White! Il988t 
Fukugita & Turned 1 99lt iKochaneU 1 9951 1 1 996at iMapz et all 
19971 IKeeton & Madaul 1200 ll ISarbu. Rusin. & Mai 1200 It 
Takah ashi & Chibal ^"ool iLi & Ostrikeil 120021 iDavis et alJ 
20031 iHuterer & Mal200llKuhlen et alJ2004tlMitcheU et alJ 
2004). Conventional wisdom holds that the statistical effects 
of ellipticity and shear are confined to the relative numbers of 
double and quadruple lenses, and that symmetric lenses are 
adequate for applications such as deriving cosmological con- 
straints. 
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To our knowledge, this conventional wisdom is based on 
a few studies in which the analysis of ellipticity and shear 
was su bordinate to pra ct ical applica t ions of l ens statistics. 
King & Browiie] il99&. iKochanekl Jl996bl) . IKeeton et all 
( 1997), and Rusin & Tegmark (2001) all computed the rela- 
tive abundances of different image configurations as a func- 
tion of ellipticity and/or shear, for various assumptions about 
the luminosity function of background sources. Along the 
way, they necessarily computed the effects of ellipticity and 
shear on the lensing cross secti on and magn ification bias, but 
did not explicitly discuss them. Chae (2003) included elliptic- 
ity in lens statistics constraints on cosmological parameters, 
but the effects were built into his statistical machinery and not 
presented on their own. We believe there is pedagogical value 
in isolating the statistical effects of ellipticity and shear and 
studying them in detail. It would be useful to lay out exactly 
how ellipticity and shear affect the lensing optical depth, and 
how that may (or may not) lead to biases in cosmological con- 
straints. Furthermore, at least two other issues deserve to be 
studied as well: the effects of ellipticity and shear on the dis- 
tribution of lens image separations; and the dependence of the 
various statistical effects on the luminosity function (LF) of 
the background sources. We will show that, while not wrong, 
the conventional wisdom is somewhat limited and there are 
effects of ellipticity and shear on lens statistics that are subtle 
but interesting. 

We focus on lensing by galaxies, by which we mean sys- 
tems with a single dominant mass component that can be ap- 
proximated as an isothermal ellipsoid. The isothermal pro- 
file describes early-type galaxies remarkably well on the 5- 
10 kpc scales r elevan t for strong lensing ( e .g.JRix et all 19971 
Gerha rd et ail 120011 iRusin & Mai l200U lTreu & Koo pmans 
120021 iKoooman s et alJl20 03: Rusin et al1 l2003l) . Lens statis- 
tics are rather different for groups and clusters of galaxies 
modeled with GNFW profiles, and tha t para llel case has re- 
cently been studied by Qguri & Keetonl ll2004l) . 

2. METHODOLOGY 
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2.1. General theory 

The probability for a source at redshift z s to be lensed can 
be written as 

t(z s )=— I dV f da -p- fdep e (e) j (f 2 7P 7 (7,</> 7 ) 



du 



$sic(Vm) 



mult 3>src(£) 
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The first integral is over the volume of the universe out to the 
source. The second integral is over the population of galaxies 
that can act as deflectors. For isothermal lenses the impor- 
tant physical parameter is the velocity dispersion, so the most 
useful description of the galaxy population is the velocity dis- 
persion distribution function (dn /da)da, or the number den- 
sity of galaxies with velo city dispersion between a and a+da 
(see Mitch eTTet al.l20 04). The third integral is over an appro- 
priate distribution p e for the internal shape (ellipticity) of the 
lens galaxy. (Without loss of generality, we can work in co- 
ordinates aligned with the major axis of the galaxy so we do 
not need to consider the galaxy position angle.) The fourth 
integral is over an appropriate distribution p 1 for the external 
tidal shear caused by objects near the lens galaxy; this inte- 
gral is two-dimensional because shear has both an amplitude 
(7) and a direction (0 7 ). Finally, the fifth integral is over the 
angular position u of the source in the source plane, and is 
limited to the multiply-imaged region. In the last integrand, 
fi is the lensing magnification, ^ SIC (L) is the cumulative num- 
ber density of sources brighter than luminosity L in the sur- 
vey, and the factor $ src (L//i)/$ slc (L) accounts for the "mag- 
nification bias" that produces an excess of faint sources in a 
flux- limited survey due to lensing magnification jTurner et al.l 
1984). (The role of the l imiting flux or limiting luminosity 
will be discussed in § !2.3l below.) The differential probability 
for having a lens with image separation Ad can be computed 
by inserting a Dirac <5-function in eq. Q to select the parame- 
ter combinations that give separation Ad. In other words, we 
can think of the (normalized) image separation distribution as 
p(A9) = r" 1 dr/dAO. 

The lensing cross section A and the magnification bias fac- 
tor B are often computed separately (see Chae 2003 for the 
most recent example). According to eq. Q, however, the two 
quantities are closely linked. We prefer to keep them together 
and compute the product 
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which we call the "biased cross section." The biased cross 
section depends on both the lens model parameters and the 
source LF. 

A convenient feature of isothermal lenses is that the physi- 
cal scale decouples from the lensing properties. All of the de- 
pendence on Zs, Zi, and a is contained in the (angular) Einstein 
radius 9e, so when we work in units of 9e nothing depends 
explicitly on these parameters. As a result, the dimensionless 
biased cross section BA = BA / 9\ depends only on the elliptic- 
ity and shear (and implicitly on the source LF). We can then 
rewrite eq. Q as 



de p e (e) / d 7 p y (j,(f>y) BA(e,i) > , (3) 



where the second factor depends only on the ellipticity and 
shear distributions, while the first factor depends only on the 
source redshift and the deflector population. If we only want 
the change in the optical depth produced by ellipticity and 
shear, then we can simply write 



= / de p e (e) / d 7J> 7 (7,<£ 7 ) 



£A(e,7> 7 ) 
BA n 



(4) 



where BAq and to are the biased cross section and the optical 
depth for the spherical case. 

Working in dimensionless units also simplifies the study 
of image separations. Even if all galaxies are SIS, the dis- 
tribution of image separations will be fairly broad because 
ther e is a range of lens galaxy masses and redshifts (see, 
e.g.. lKochanekl[l993al) . However, the dimensionless separa- 
tion A§ = A0/6 E is always AO = 2 for an SIS lens, so the 
dimensionless image separation distribution p(A8) is just a 6- 
function when the ellipticity and shear are zero. This means 
that studying p(A9) is the simplest way to identify changes 
to the image separation distribution caused by ellipticity and 
shear. For fixed ellipticity and shear, the distribution can be 
formally written as 



p(A9\e, 7 ,<t>J-- 



du 
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where A$(u; e, 7, 4>~) is the dimensionless image separation 
produced for a source at position u by a lens with the speci- 
fied ellipticity and shear. The full image separation distribu- 
tion can then be found by integrating over appropriate distri- 
butions of ellipticity and shear. Note that we do not actually 
need to integrate over the masses and redshifts of the deflector 
population in order to compute changes to the optical depth 
and the image separation distribution. 

2.2. The isothermal ellipsoid with shear 

We first discuss isothermal ellipsoids without an external 
shear, and then discuss properties of shear at the end of this 
subsection. 

The projected surface mass density for an isothermal ellip- 
soid, written in polar coordinates (r, tp) and expressed in units 
of the critical density for lensing, is 



K(r, <j>) = ■ 



1+q 2 



2r L(l+<T)-(l-<f)cos2(/) 



1/2 



(6) 



where q < 1 is the axis ratio, and the ellipticity is 
e = 1 - q. (Recall that we are working in coordinates 
aligned with the major axis of the galaxy.) The ra- 
dius r and the parameter b both have the dimensions of 
length, and may be expressed as physical lengths (e.g., 
kiloparsecs) or angles on the sky (radians or arcseconds); 
we work in angular units. The le nsing properties of 
an iso th ermal ellipsoid are given bv iKassiola & Kovneri 
\ 1993h. iKormann. Schneider. & B artelmannl i ll 9941) . and 
Keeton & Koc hanekll|1998l) . 

The Einstein radius sets the lensing scale, so it is useful to 
determine its value. Consider the deflection anM produced 
by the monopole moment of the lens galaxy, 



ao(r) ■■ 



1 [dr> r d( f>r' K (r'^ = ^l, (7) 
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where M cy \(r) is the projected mass in a cylinder of radius r. 
The Einstein radius is defined by oo(&e) = Oe- This defini- 
tion reduces to the standard Einstein radius in the spherical 
case, and it is the quantity that seems t o be most relev ant in 
models of nonspherical lenses (e.g.. ICohn et alJ [2001 ). For 
the isothermal ellipsoid, we find 



b 



1 



[2(1 +q- 2 )] 1/2 K(l + q 



(8) 



where K(x) is the elliptic integral of the first kind. For refer- 
ence we note that this function can be approximated by 



e E (e) 



= exp [(0.89e) 3 



(9) 



which is accurate to <1% fore < 0.53 and to <4% fore < 0.9. 
In practice, however, we use the exact result. 

We must specify how to normalize the model, or how to 
choose the parameter b. For a spherical galaxy, b simply 
equals the Einstein radius and is related to the velocity dis- 
persion by 



\ c 



D„ 



(10) 



where D os and D/ s are angular diameter distances from the 
observer to the source and from the lens to the source. For a 
nonspherical galaxy the situation is less straightforward. If 
we seek a dynamical normalization in terms of a measur- 
able stellar velocity dispersion, then we must worry about 
complications invol ving the halo shape and proj ection effects 
JKeeton et al.lli997t iKeeton & Kochanekl 119981 IChael 2003). 
Consider the dynam ical normalization shown in Figure[2(fol- 
lowing Chae 2003). At a typical ellipticity e w 0.3, b could 
rise by 7% (compared to the spherical value) if all halos are 
oblate, or fall by 7% if all halos are prolate. Some dissipation- 
less numerical simulations have predicted roughly compara- 
ble numbers of oblate and prolate halos (Dubinski & Carlberg 
1991; Jing & Suto 2002), which would yield a b value less 
than 1% higher than the spherical value (for e = 0.3). How- 
ever, the shape distribution is likely to be affected by hydrody- 
namics (e.g.. lKazantzidis et al.ll2004l) . so it is not understood 
in detail (in simulations, let alone in reality). In other words, 
the dynamical normalization appears to be small but uncer- 
tain, and impossible to compute precisely. 

An alternate approach is to fix the Einstein radius to be 
independent of ellipticity (and shear; see below). This 
seems reasonable, because the Einstein radii of observed 
lenses can generally be determined in a model-independent 
way to a few percent accuracy (e.g., ICohn et ail 1200 It 
Munoz, Kocha nek. & Keetonl2 001). and because it keeps the 
mass properties (the aperture mass) independent of ellipticity 
and shear. This normalization is also shown in Figure^ and it 
is the one we adopt. However, it is important to keep in mind 
that there is an irreducible uncertainty of a few percent in our 
analysis associated with the normalization. 

Objects in the vicinity of the lens galaxy create tidal forces 
that affect the lens potential. The contribution is often mod- 
eled as an external shear whose contribution to the lens poten- 
tial is 

(11) 



$shear(r, <P) = ~ ^ COS 2(0- <f>~) , 



where (r, cf>) are polar coordinates, 7 is the dimensionless shear 
amplitude, and 7 is the direction of the shear. As an example, 
consider the shear produced by an isothermal sphere galaxy 
with Einstein radius bo that lies at polar coordinates (ro,</>o) 
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FIG. 1 . — Change in the isothermal ellipsoid b parameter as a function of 
ellipticity. The two dashed curves show the dynamical normalization if all ha- 
los are assumed to be oblate or prolate. The dot-dashed curve shows the case 
when half of the halos are assumed to be oblate and half prolate. The solid 
curve shows the result when the Einstein radius is fixed to be independent of 
ellipticity, which is what we assume in this paper. 



relative to the lens galaxy; the shear amplitude would be 7 = 
Z?o/(2ro) and the shear direction would be </> 7 = <f>Q. External 
shear does not contribute to the local surface mass density, 
so it does not affect the monopole deflection or the Einstein 
radius. 

2.3. Source luminosity functions 

The number density of sources with luminosity be- 
tween L and L + dL is given by the luminosity function 
[c/0 slc (L) /dL]dL. The quantity of interest for lens statistics 
(see eq. 0) is the cumulative number density of sources 
brighter than L, or 



dL 1 



dL 1 . 



(12) 



We consider model LFs appropriate to both radio and optical 
surveys. 

The simplest model LF is a featureless power law, <p S ic(L) oc 
L~P . In this case the biased cross section simplifies to 



BA ■■ 



/9-1 



mult 



(13) 



where p(p) is the distribution of magnifications for lensed 
sources. Several points are worth mentioning. First, with 
f3 — > 1 + the magnification weighting factor becomes unity and 
we recover the simple lensing cross section with no magni- 
fication bias. 8 Second, because the magnification distribu- 
tion generi cally has a power law tail p(u) oc /z~ 3 for /i ^> 1 
(see Schneider, Ehlers, & Falco 1992), the integral diverges 
for f3 > 3 and the biased cross section is well defined only 
for f3 < 3. Finally, a power law is featureless so the bi- 
ased cross section does not depend on the particular flux or 
luminosity limit of a survey. A power law LF is a good 
model for radio surveys. For example, the largest existing 
lens survey is the JVAS/CLASS survey of flat-spectrum radio 

8 We write (3 — » 1 + , meaning that j3 app roac hes unity from above, because 
for /9 < 1 the cumulative LF integral (eq. G3) formally diverges. Neverthe- 
less, the biased cross section integral (eq. 1131 ) remains well defined. 
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sources ( Mvers et all2003tlBrowne et all2003l) . and it has an 
LF that is well d escribed by a pow er law with /3 w 2.1 (see 
Rusin & Tegm arkfeOOH IChael2003T) . 

Future lens samples are likely to be dominated by optically- 
selected quasar lens es found in deep wide-field imaging sur- 
veys (e.g.. iKuhlen et all 120041) . While accurate determina- 
tion of the quasar LF is a long-standing problem, recent 
evidence favors the double pow er law form proposed by 
iBovle. Shanks. & Peterson! i ll 9881) . 

d<j)(L,z) ^* dL_ 

dL [L/U(z)W'+[L/U(zW> U(z) ' 

where the break luminosity L* evolves with redshift as 
jMadau et al.ll999t) 

L*(z) = L,(0) (1 +zT s "' j~ Tz , (15) 

where the quasar spectral energy distribution is assumed to be 
a power law, /„ oc v~ a <- . With this LF, the biased cross sec- 
tion clearly depends on the bright and faint slopes f3h and /?/, 
and also on the limiting luminosity L cut (z)/L*(z). It depends 
on source redshift to the extent that the se quant i ties d epend 
on redshift. We adopt the model from Fan et al. (20QJJ) with 
bright-end slope /3/, = 3.43 at z < 3 and (3h = 2.58 at z > 3, 
and fa int-end slope /?/ = 1 .64 at all redshifts (IWvithe & Loebl 
2002). 

If we wanted to compute statistics for real quasar lens sur- 
veys, we would need to adopt an appropriate limiting mag- 
nitude and compute the limiting luminosity L cut (z) / L*(z) as a 
function of redshift. This would require specifying the pass- 
band, computing /f-corrections, and other details that would 
muddy the waters. Since the goal is conceptual understanding 
of the effects, we believe that it is simpler and more instruc- 
tive to work with a luminosity cut L cut /L*. In this case, we do 
not need to specify the parameters c/>*, L*(0), z*, C> an d £■ 

2.4. Numerical techniques 

We compute the integrals in eqs. (|2j-(|3 using Monte Carlo 
techniques. First, for fixed ellipticity and shear we place 10 5 - 
10 6 random sources in the source plane, in the smallest circle 
enclosing the causti cs. We solve the lens equation using the 
gravlens software ( Keeton 2001) to determine the number of 
images and their positions and magnifications. We define the 
image separation to be the maximal separation between any 
two images in the system, Ad = max \8i-0j\; this is a conve- 
nient, observable, and well-defined quantity that is indepen- 
dent of the number of images. 

We separate the lenses into three standard classes based on 
the image multiplicity: "doubles" have two bright images, one 
with positive parity and one negative, plus a faint central im- 
age that is rarely observed; "quads" have four bright images, 
two positive and two negative parity, plus a faint central image 
that is rarely observed; and "naked cusps" have three bright 
images, either two positive parity and one negative or vice 
versa. We use the classifications directly only when studying 
the quadruple-to- double ratio (§|5}. The classification offers 
a fringe benefit: we can identify numerical errors as systems 
that do not fit into any of the classes (because, for example, 
the software failed to find one of the images). We estimate 
that the numerical failure rate is < 10~ 4 . 

Next, where appropriate we integrate over distributions for 
ellipticity and shear. For the ellipticity, we adopt the distribu- 
tion of ellipticities measured for 379 early-type galaxies in 1 1 




FIG. 2. — Magnification distributions for spherical deflectors without shear 
and with shear 7 = 0.2. The curves are normalized so that the area under each 
curve is the corresponding cross section in units of d\. For the shear case we 
show the distributions for doubles and quadruples separately. Note that the 
minimum magnification for doubles is /i2,min = 2/[(l + 37XI-7)], while for 
quadruples it is /M.min - 2/[7(l~7 2 )] IFinch et alJ20Q2l) : so /i2.min = 1 -56 and 
M4,min = 10. 4 for the case 7 = 0.2 shown here. The distributions asymptote to 
A(/i) oc /i -3 at high magnifications 1 Schneider et al. 1992). 



nearby clusters by J0rg ensen et alJ (I1995I) . The distribution 
has mean (e) = 0.31 and dispersion a e = 0.18, and there are 
no galaxies with e > 0.8. Although the measured ellipticities 
describe the luminosity while what we need for lensing is the 
ellipticity of the mass distribution, this is probably the best 
we can do at the moment. In any case, it seems reasonable 
to think that the ellipticity distribu tions for the lig ht and the 
mass may be similar (see iR usin & Tegmark 120011) . For the 
shear, Holder & Schechter (2003) estimate that the distribu- 
tion of shear amplitudes derived from simulations of galaxy 
formation can be described as a lognormal distribution with 
median 7 = 0.05 and dispersion c 7 = 0.2 dex; this distribution 
is also broadly consistent with the shears required to fit ob- 
served lenses. As a rule of thumb, a shear 7 ~ 0. 1 is common 
for lenses in poor groups of galaxies, and the shear can reach 
7 ~ 0.3 for lenses in rich clusters (e.g., j Keeton et alJ|1997t 
iKundic et alJl997albtlFischer et alll998HKneib et all2000l) . 
We assume random shear orientations. 

3. THE OPTICAL DEPTH 

Before determining the effects of ellipticity and shear on 
the lensing optical depth, it is instructive to consider first how 
they affect the source plane. There is only a small change in 
the lensing cross section. In fact, for isothermal galaxies shear 
has no effect on the radial caustic and hence on the cross sec- 
tion. 9 Ellipticity (or any other internal angular structure) in 
isothermal galaxies changes the caustics in such a way as to 
reduce the cross section, as explained in the Appendix. The 
main effect of increasing ellipticity or shear is to lengthen the 
tangential caustic, which enlarges the phase space for large 
magnifications and raises the tail of the magnification dis- 
tribution, as illustrated in Figure [2] In particular, we see a 
sharp increase in the cross section for producing magnifica- 

9 Shear can affect the cross section only in the rare case that the tangential 
causti c pierces the radial caustic to form naked cusps (e.g., Schne ider et all 
1992). For SIS+shear models, this happens only when the shear is large, 
7 > 1/3. In this case, there is some (small) multiply-imaged region outside 
the radial caustic. 
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tions larger than the minimum magnification for a quadruple 
lens (see caption). 

We now examine the dependence of the biased cross sec- 
tion BA on ellipticity (Figure[3j and shear (Figure|4}. When 
the LF is a power law with (3 — > 1 + there is no magnifica- 
tion bias, and Figure [3] illustrates how ellipticity reduces the 
cross section. Even with magnification bias, ellipticities up to 
e ~ 0.5 do not affect the biased cross section by more than 
10% unless the source LF is very steep (e.g., the very bright- 
est quasars, L cut /L* > 100 when /?/, = 3.43). Figure |4] shows 
that shear causes a stronger increase in the biased cross sec- 
tion, but we must remember that realistic shears are 7 < 0.1 
and only lenses in clusters experience large shears of 7 ~ 0.2- 
0.3. Thus, the typical change in the biased cross section due to 
shear is again no more than 10% unless the LF is very steep. 

To compute changes in the full optical depth we integrate 
the biased optical depth over appropriate distributions of el- 
lipticity and shear (see eq. @). The results are shown in Fig- 
ure Without magnification bias {(3 — > 1 + ), ellipticity and 
shear reduce the optical depth very slightly (by ~0.6%, al- 
though the statistical uncertainty from our Monte Carlo cal- 
culations is ~0.3%) With magnification bias and a shallow 
source LF, ellipticity and shear can reduce the optical depth 
by up to ^2.5% relative to the spherical case. For power law 
LFs, only when j3 > 2.2 is there an increase in the optical 
depth, and we must have (3 > 2.5 (J3 > 2.6) in order for the 
increase to be more than 5% (10%). For the quasar LF, the 
increase exceeds 5% only if the bright end is steep (J3/, = 3.43) 
and the survey is limited to bright quasars (L cut /L* > 10). 

In practice, these results mean that ellipticity and shear are 
important for the optical depth only in surveys that are re- 
stricted to the brightest quasars. They are not very significant 
for the sorts of deep optical surveys now underway that probe 
well beyond the break in the quasar LF. 

4. THE IMAGE SEPARATION DISTRIBUTION 

We now turn to the distribution of lens image separations 
and how it is affected by ellipticity and shear. First, we re- 
call several basic facts. Even in spherical models the distribu- 
tion of dimensioned image separations will have some natural 
spread because of the range of galaxy masses and redshifts; 
but the distribution of dimensionless separations AO is a 6- 
function at AO = 2. To highlight changes in the separation 
distribution, it is therefore useful to focus on the distribution 
of Ad. Also, as discussed in § 12.41 we define the separation to 
be the maximal distance between any pair of images. 

Figure[6]shows the distribution of AO for several values of 
ellipticity (upper panel) and shear (lower panel). The distri- 
bution has an interesting shape that peaks at the ends and is 
low in the middle. It has a sharp cutoff at the high end, while 
at the low end it has a sharp drop followed by a small tail to 
lower values. The peaks correspond to sources near the minor 
and major axes of the lens potential. 

As the ellipticity or shear increases, the distribution of AO 
broadens and its mean shifts. To quantify these effects, we 
compute the mean separation (AO) and the spread a A g = 

((AO 2 ) - (AO) 1 ) 1 ! 2 , and plot them as a function of elliptic- 
ity or shear in Figure The increase in the mean and scatter 
are small for all ellipticities, and are both <20% for all but 
the strongest shears (7 > 0.3) felt by lenses in cluster envi- 
ronments. Nevertheless, it is interesting that shear produces a 
net bias toward larger image separations. 
Finally, by averaging over the ellipticity and shear distribu- 



tions we obtain the net image separation distribution shown in 
Figure [8] The averaging has smoothed out the sharp features 
seen in Figure [6] when the ellipticity and shear were fixed. 
The net distribution is nearly Gaussian, with mean AO = 2.01 
and scatter er A( j = 0.18 for a power law LF with (3 = 2.1, or 

AO « 2.01 and a A g « 0.19 for various cases of the quasar LF. 
In other words, ellipticity and shear basically leave the mean 
image separation unchanged but create an additional scatter 
of 10%, and these results are insensitive to the source LF. 

5. QUADRUPLE-TO-DOUBLE RATIO 

We next consider how ellipticity and shear affect the num- 
ber of lenses with different image configurations. While an 
SIS lens always produces two images, increasing ellipticity or 
shear leads to increasing probability for configurations with 
four images. Furthermore, large ellipticities (e > 0.606) or 
shears (7 > 1/3) can lead t o "naked cusp" co nfigurations 
with three bright images (e.g. JKeeton et alll997l) . Nearly all 
known lenses with point-like images are doubles or quadru- 
ples; among ~80 known lenses the re is only one can didate 
naked cusp lens (APM 08279+5255 : lLewis et alJ2002l) . 

Figure [9] shows that the quadruple to double ratio rises 
monotonically with ellipticity or shear. For the CLASS LF, 
the ratio is ^20% for typical ellipticities e ~ 0.3 or shears 
7 ~ 0.1. Our results agree well wit h previous analyses (e.g., 
iRusin & TegmarU200lllFinch et all2002l) . 

We can now estimate the expected number of quadruples 
(and cusp triples) by averaging over our fiducial distribution 
of ellipticity and shear. The results are shown in Figure [K)] 
for both the power law LF and the quasar LF. For the CLASS 
LF (f3 = 2.1), the net quadruple-to-double ratio is about 0.35, 
while the triple-to-double ratio is only about 0.01. The num- 
ber of quadruple vs. double systems in the CLASS statistical 
sample is 5 vs. 7; the ratio is t wice as large as our pred ic- 
tion. We therefore agree with Rusin & Tegmark (2001) in 
concluding that ellipticity and shear alone cannot easily ex- 
plain the high number of quadruples. Additional effects are 
required, which are probably related to lens galaxy environ- 
ments. Shear is only a low-order approximation to the lensing 
effects of objects near the lens galaxy. Recent studies have 
sho wn that including higher-order effects from satellite galax- 
ies JCohn & Kochanekl 2003 ) or extended groups of galaxies 
(Keeton & Zabludoff 2004) around the lens can significantly 
boost the quadruple-to-double ratio. 

Figure flOl shows that surveys targeting lensed quasars are 
expected to have a low quadruple-to-double ratio unless the 
bright end of the LF is steep (/?/, = 3.43) and the survey is lim- 
ited to bright quasars (L cut /L* > 10). This prediction could, 
of course, be an underestimate because we have neglected the 
higher order effects from lens environments. 

6. EFFECTS ON COSMOLOGICAL CONSTRAINTS 

While the changes in the optical depth and image separa- 
tion distribution caused by ellipticity and shear seem mild, it 
is important to quantify how they affect one of the main appli- 
cations of lens statistics: constraints on cosmological param- 
eters. One approach would be to modify the analyses of real 
lens samples to include the full effects of ellipticity and shear 
(building upon the analysis of Chae 2003). Such an approach, 
however, would be limited by Poisson uncertainties in current 
lens samples (e.g., CLASS has just 13 lenses), by systematic 
uncertainties where models may or may not b e correct (e.g., 
evolu tion in the lens galaxy population; see Mitchell et al. 
2004), and by systematic effects that are known to be present 
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FIG. 3. — Enhancement in the biased cross section caused by ellipticity. The shear is set to zero, (a, left) Results for different power law source LFs. Recall 
that with /3 — ► 1 + there is no magnification bias, and that the CLASS radio survey has fa 2.1. (b, right) Results for the model quasar LF, for different values of 
the limiting luminosity (L cllt /L*) and the bright-end slope (/3;,); the faint-end slope is fixed at /3/ = 1.64. The jaggedness in the upper curves is due to statistical 
noise, because for steep LFs the magnification bias is strong and the results are dominated by rare extreme-magnification systems. We estimate the statistical 
errors in the upper curves to be ~3%, and much smaller for the other curves. 
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FIG. 4. — Similar to Figure^] but for shear. The ellipticity is set to zero here. The statistical errors in the upper curves are ~10%, and much smaller for the 
other curves. 
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FIG. 6. — Histograms of the dimensionless image separation AS = A9/9e, 
for two values of ellipticity (top panel) and shear (bottom). For a spheri- 
cal lens the distribution is a (5-function at A8 = 2, indicated by the vertical 
line. The histograms contain all image multiplicities (i.e., both doubles and 
quadruples). The results are shown for the CLASS LF (a power law with 
/3 = 2.1), but they are not very sensitive to this choice. 
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FIG. 7. — Mean and spread of the image separation distribution as a func- 
tion of ellipticity or shear. The image separations are in units of Be- As in 
Figurelrjl the results are shown for the CLASS LF but are not very sensitive 
to this choice. 



in the data but have not yet been studied (e.g., having multiple 
lens galaxies). We believe that it is more instructive to create 
mock lens surveys that mimic CLASS but allow us to isolate 
the effects of ellipticity and shear. Specifically, we create sur- 
veys that include ellipticity and shear, and then analyze them 
using standard spherical models in order to uncover biases 
that result from neglecting ellipticity and shear. We create 
mock surveys with 1000 lenses in order to minimize Poisson 
uncertainties. We use a Monte Carlo approach, drawing pa- 
rameter values from appropriate probability distributions (as 
indicated in eq.Q. Specifically: 

• A subset of sources in the CLASS survey has a redshift 
distribution that can be treated as a Gau ssian with mea n 
(z s ) = 1.27 and width a z = 0.95 JMarlowet aU feOOO). 
and this is usually taken as a mod el for the re dshift 
distribution of the full survey (e.g.. lChaell2003l) . The 
Gaussian is modified by the redshift dependence of the 




FIG. 8. — Net image separation distribution after averaging over ellip- 
ticity and shear. The histogram is normalized to unit area. The results are 
again shown for the CLASS LF but are not very sensitive to this choice. The 




FIG. 9. — Quadruple-to-double ratio as a function of ellipticity or shear, 
assuming the CLASS LF (a power law with /3 = 2. 1). 



optical depth t(z s ) to obtain the redshift distribution of 
the lensed sources (see Mitchell et all2004l) . 

• The lens redshift is drawn from the distribution p(zi) oc 
(Di s /D os ) 2 dV/dzi, where the factor of (D ts /D os ) 2 
comes from the factor of Q\ in the lensing cross sec- 
tion (see eq. (I10» . 

• The velocity dispersion is drawn from p(a) tx 
u dn/da, where the factor of <r 4 comes from d\ in 
the lensing cross section. We use the velocity disper- 
sion d istribution function dn /da derived by Shet h et ail 
(2003) for early-type galaxies in the Sloan Digital Sky 
Survey. For simplicity, we assume that the velocity dis- 
persion distribution does not evolve with redshift. We 
could add evolution to both the creati on and an alysis 
of the mock survey (see, e.g., Mitchell et al. 2004), but 
that would just complicate matters. 



We us e the ellipticity distribution from Jcirg ensen et alJ 
(119951) . the shear amplitude distribution from 
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FIG. 10. — Quadruple-to-double ratio, as a function of the slope (i of power law LFs (left) or the limiting luminosity L cut /L t for quasar LFs (right). The results 
are obtained by averaging over the distributions of ellipticity and shear discussed in the text. In the left panel, we also show the ratio of triple (or naked cusp) 
lenses as the dashed curve. 



iHolder & Schechteil J2003I) . and random shear di- 
rections. 

• When drawing random source positions, we use magni- 
fication bias appropriate to the CLASS survey (a power 
law with (3 = 2.1) since it is the most commonly used 
survey in current lens statistic analyses. 

Given the parameters we can compute the observables for 
each mock lens: the source and lens redshifts and the image 
separation. The other key observable is the total number of 
sources in the survey. We use the optical depth to determine 
the number of sources needed to obtain 1000 lenses, which 
is typically ^8 x 10 5 . We distribute these sources in redshift 
using the Gaussian given above. 

We then analyze the mock survey with standard maximum 
likelihood techniques. Assuming complete data — knowl- 
edge of the image separation and the lens and source redshifts 
for lens systems, and the redshift distribution of non-lensed 
sources — we use the likelihood function 



£ = 



(N pKd f 



\ prcil 



N ohs l 



n 



i 



d 2 



t(z s ,0 dzu dAdi 



(16) 



The first term represents the Poisson probability for having 
A'obs observed lenses when N pK d are predicted, while the sec- 
ond term represents the probability that the lenses have the 
observed properties (e.g., observed lens redshift zi and image 
separation A9 given the source redshift z s ). As mentioned 
above, we neglect ellipticity and shear in the likelihood anal- 
ysis because we want to understand the biases that may occur 
when spherical models are used for lens statistic analyses. We 
hold the parameters in the velocity dispersion function fixed 
at their input values si nce unc ertainties in these parameters 
have negligible effect (Mitch ell et a 1 J 12004), Thus, the only 
variables in the model are the cosmological parameters Q,m 
and which we adjust to maximize the likelihood. We use 
input values of VLm = 0.3 and ft\ = 0.7, and study how much 
the recovered values differ. As mentioned above, using sur- 
veys with 1000 lenses should mitigate Poisson uncertainties, 
but we always produce and analyze 10 independent surveys to 
verify that the statistical noise in our results is negligible. 

It is useful to begin by examining two toy models that fo- 
cus on how changes in the optical depth or image separation 
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FIG. 11. — Biases in constraints on cosmological parameters from 
analyses of lens statistics. We show the errors Af2« = f^™ d _ an ^ 
AOa = f2'? 0<1 fiA Ue that result from using simple models with spherical lenses 
that neglect shifts in the optical depth (crosses), image separations (triangles), 
or both. (See text for details.) The statistical uncertainties are smaller than 
the size of the points. 



distribution can affect cosmological constraints. In the first 
case, we imagine using spherical lens models but manually 
adjusting the optical depth. This is equivalent to changing the 
total number of deflectors. In practice, it means adjusting the 
number of sources in our mock survey (since we fix the num- 
ber of lenses). The crosses in Figure ^2 show the errors in 
the recovered cosmological parameters if the difference be- 
tween the actual (input) optical depth and the spherical model 
is (-20,-15,-10,-5,0,5,10,15,20)%. We see that simply 
changing the optical depth moves the cosmological parame- 
ters mainly along the line corresponding to flat cosmologies, 
and the shift is fairly small: A51a = 0.03 if the real optical 
depth is 10% larger than predicted by the spherical model. 
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In the second case, we again start with spherical models 
but manually adjust the image separations. This is equiva- 
lent to shifting the velocity dispersion distribution to higher 
or lower values, and then adjusting the number of galax- 
ies to keep the optical depth fixed. The triangles in Fig- 
ure!^ show the results of shifting the image separations by 
(-20,-15,-10,-5,0,5, 10, 15,20)%. There is a large shift in 
the recovered cosmological parameters, and it is almost or- 
thogonal to the line of flat cosmologies. For example, if the 
real image separations are 10% larger than predicted by spher- 
ical models, then there will be errors of A£Im = 0.12 and 
AOa = 0.21 in the parameters recovered by spherical mod- 
els. These two cases are just toy examples, but they illustrate 
the important principle that even small errors in the model im- 
age separations can have a significant effect on cosmological 
constraints (even if small errors in the optical depth do not). 

Finally, we consider the case where we use the full effects 
of ellipticity and shear on the mock survey. Essentially, this 
amounts to using the corrections to the optical depth from 
Figure [5] and to the image separation distribution from Fig- 
ure [8] The circle in Figure ITT1 shows that neglecting ellip- 
ticity and shear in the likelihood analysis causes errors of 
AO M = 0.00 ± 0.01 and Af2 A = -0.02 ± 0.01, where the er- 
rorbars represent the statistical uncertainties in our calcula- 
tions. (We have achieved small Poisson uncertainties but 
not eliminated them altogether.) That is the case if we al- 
low Hm and fl\ to vary independently. If we restrict atten- 
tion to flat cosmologies (Hm + ^a = 1) then the bias is just 
Afijf = — AOa = 0.01 (with negligible errorbars). This result 
is consistent with our conclusions from the previous sections 
that ellipticity and shear have little effect on the optical depth 
and mean image separation. It is nonetheless valuable to have 
a careful validation of the conventional wisdom that ellipticity 
and shear do not significantly affect cosmological constraints 
derived from lens statistics. 

7. CONCLUSIONS 

The effects of ellipticity and shear on strong lensing statis- 
tics have been swept unde r the rug in most analyses to date 
(a valiant exception being Chae 2003). The reason for this 
is twofold: (1) models with nonspherical deflectors intro- 
duce new, and sometimes poorly constrained, parameters and 
greatly complicate calculations; and (2) conventional wisdom 
suggested that realistic ellipticities and shears had little effect 
on anything but the image multiplicities. We have stepped be- 
yond the state of blissful ignorance to present a general anal- 
ysis of how ellipticity and shear enter into lens statistics. 

The effects depend strongly on magnification bias, which 
in turn depends on the luminosity function of sources in a 
lens survey. If the LF is a power law oc L~^, as in radio sur- 
veys like CLASS (with (3 rj 2.1), ellipticity and shear gen- 
erally decrease or increase the lensing optical depth by only 
a few percent. The increase is more than 5% only if the LF 
is steep (J3 > 2.5). For optical quasar surveys, if the limit- 
ing luminosity is below the break in the quasar LF then ellip- 
ticity and shear decrease the optical depth by a few percent. 
There is a noticeable (>5%) increase only for surveys limited 
to the brightest quasars (L cut /L* > 10, if the bright end slope is 
(3h = 3.43). Since ongoing and planned optical surveys are ex- 
pected to reach to the faint end of the QSO LF (L cut /L* < 1), 



we do not expect ellipticity and shear to have a large effect on 
the predicted number of lenses in future lens surveys. 

Ellipticity and shear do not shift the mean of the distribution 
of lens image separations, but they do introduce an additional 
scatter of ^10%. They naturally affect the relative numbers 
of double, quadruple, and triple lenses, but they cannot easily 
explain the high observed quadruple-to-double ratio. Ellip- 
ticity has little effect on predictions for elusive central lensed 
images, although it does lead to a segregation that quadruple 
lenses are generall y expected to have fainter central images 
than double lenses (Keeton 20031). 

Since ellipticity and shear produce only small changes in 
the lensing optical depth and image separation distribution, 
they are not very important in lensing constraints on cosmo- 
logical parameters. Neglecting them leads to biases in £Im 
and £1a of <0.02. Moreover, hydrodynamical N-body simu- 
lations tend to find systems that are more spherical than thos e 
in dissipationless simulations (e.g.. iKazantzidis et alJl2004l) . 
Therefore, the ellipticity effects on lensing statistics found in 
this paper, while already small, could even be an overestimate. 

We conclude that for lens statistics problems other than im- 
age multiplicities, ellipticity and shear have surprisingly little 
effect. Unless percent-level precision is needed, or a survey 
with a particularly steep LF is being considered, ellipticity 
and shear can probably be ignored. Their effects will be- 
come more important as lens samples grow into the hundreds 
or thousands and s tatistical uncertainties plummet (see, e.g., 
iKuhlen et all2004T) . At that time it will be important to know 
the distributions of ellipticity and shear, and also to resolve 
questions about how to normalize the lens models (see § !2.2t . 

There are systematics besides ellipticity and shear that 
may affect strong lens statistics. They include merge rs and 
evolutio n in the deflecto r population (e .g., |Rjx_£LaL 1994; 
| Mao & Kochanekl 1 19941 iKeetonl [2002L JOfek et all 120031: 
IChae & Mad 120031 iMitchell et alJ 12004 . halo triaxiality 
(e.g., Oguri & Keeton 2004) or ot h er complex internal 
structure (e.g.. IMolleretalJ 120031 lOuadrietalJ [ 2003). 

comp o und lens galax ies (e.g. . iKo chanek & Apostolakis 

Il988l iMoller & Blair] l200lt ICohn & Kochanekl l2003l)T 
and lens galaxy environ ments (e.g., IMolleretalJ 2002; 
Keeton & Zabludoff 2004). In order to bring lens statistics 
into the realm of precision cosmology, each of these factors 
must be addressed carefully. We have taken one step in that 
direction by studying ellipticity and shear, finding that their 
effects are relatively small and in principle easy to take into 
account. 
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APPENDIX 
WHY IS A/A < 1? 

In this appendix we derive the cross section for a generalized isothermal lens to explain the result from §[3] that ellipticity 
reduces the cross section. The lens potential for a generalized isothermal model has the form <f>i so (r, <f>) = r /(</>) where f(<p) is an 
arbitrary function specifying the angular shape. Consider expanding the potential in angular multipoles, 



*-Ur,<i>) = E r 1 



OO 

■n 

m=l 



a„, cos(m0) + b m sm(m<fi) 



where Oe is the Einstein radius (as defined in § 12.2k The corresponding mass distribution then has the form 



KisoO, </>) = ( 1 e m cos[ot(0- cj) m )] ) , 



2r 

\ m=L 

where the amplitude e,„ and direction <p m of the mass multipole are given by 

e ln = (tn 2 -l)^Jaf n + bl , 

, l-i b m 
(p m = — tan — . 

m a m 



(Al) 

(A2) 

(A3) 
(A4) 



In other words, we can think of this model as a multipole expansion of the surface mass density. 

The radi al caustic — properly term ed a pseudo-caustic since a singular isothermal lens does not formally have a radial critical 
curve (see Evans & Wilkinson 1998) — can then be written in parametric form as: 



/ OO 

«caus(A) =—9e I cosA + ^^ |(a m cos mX + b m sin ra A) cos A + m(a m sin mX — b m cos mX) sin A 

\ m=l 

(oo 
sin A + j^(a m cos m X + b m sin mX) sin A — m(a m sin m X — b m cos m A) cos A 
m=l 



(A5) 



(A6) 



Although this form appears complicated, if we collect the two coordinates w caus and v caus into a vector u caus then we can easily 
evaluate the area inside the radial caustic: 



A = 



2tt 



Ucaus (A) x 



dU c 



dX 



dX=TT0l 



1 OO 



m=\ 



(A7) 



This is the lensing cross section (provided there are no naked cusps). The summand, and hence the sum, is manifestly nonnegative, 
so the cross section for any nonspherical model is A < Aq = ir6 E . This result is illustrated in Figure lATSl for different multipole 
terms. It is clear that asphericity deforms the caustics in such a way that the cross section is smaller than for the spherical case. 
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